function  centroids_subpixel = detect_beads_LS_2(LS)

    % 预处理：高斯滤波平滑体数据，去噪
    volFiltered = imgaussfilt3(LS, 2); % 3D高斯滤波
    
    % 二值化处理
    threshold = graythresh(volFiltered(:)); % 自动阈值
    bwVol = imbinarize(volFiltered, threshold);
    
    % 去除小噪声（3D连通区域分析）
    bwVol = bwareaopen(bwVol, 50); % 删除小于50体素的连通区域
    
    % 提取3D中心坐标
    stats = regionprops3(bwVol, 'Centroid'); 
    
    centroids_subpixel = table2array(stats);
end
     